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Background:  While  metamaterials  enable  unprecedented  control  over  propagation  of  light  with 
applications  such  as  lenses  beating  the  diffraction  limit  for  hyperfine  imaging  and  lithography, 
fabrication  of  metamaterials  for  the  optical  regime  is  extremely  difficult,  as  highly  structured, 
repetitive  metallic  features  at  the  nanoscale  are  required.  The  main  barrier  to  experimental 
realization  of  optical  metamaterials  is  the  great  difficulty  of  fabricating  on  this  scale.  Approaches 
divide  into  two  main  classes.  Bottom-up  approaches  include  molecular  self-assembly  and  self¬ 
ordering.  Whilst  having  the  advantage  of  self-assembly,  the  structures  are,  as  a  result,  typically  quite 
limited  as  it  is  difficult  to  introduce  spatial  variations  required  for  most  devices.  Top-down 
approaches  such  as  e-beam  or  focused  ion  beam  lithography  can  readily  produce  complex 
structures,  however  they  have  to  individually  create  every 
part  of  the  structure,  which  is  a  time  consuming  use  of 

extremely  expensive  equipment,  and  limited  to  shallow,  (“1  Metal-dielectric  ^ 
mostly  2D,  structures.  No  existing  technology  is  suitable  preform 

to  produce  materials  structured  at  nanoscopic  scales  in 
bulk. 
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Vision:  This  project's  vision  is  to  create  a  novel  and 
practical  nanofabrication  technique  for  metamaterials  for 
frequencies  ranging  from  the  far  infrared  to  visible  light.  Figure  l:  Princinple  of  the  drawing  technique 
This  technique  will  make  inexpensive  volume  production 

of  metamaterials  possible  for  the  first  time,  It  will  do  this  by  extending  drawing,  the  technology  used 
for  making  optical  fiber,  to  preforms  containing  both  metals  and  dielectrics. 


Results  Summary:  Concentrating  on  the  more  accessible  THz  spectrum,  we  have  demonstrated  that 
the  technique  is  viable,  by  fabricating  and  characterizing  both  electric  [1]  and  magnetic 
metamaterials,  including  materials  with  measured  negative  permittivity  [2,3].  The  electric 
metamaterials  are  made  by  drawing  polymer  preforms  containing  indium  wires  (Figure  1),  and  we 
were  able  to  experimentally  demonstrate  that  this  arrangement  behaves  as  plasmonic  materials 
with  plasma  frequency  in  the  Thz  spectrum  [1].  Magnetic  materials  are  more  challenging  to  fabricate 
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as  they  require  non-circular  metallic  wires,  which  we  were  able  to  produce  using  two  separate 
method.  The  first  method  on  which  we  worked  in  collaboration  with  our  colleague  M.  Large,  was  to 
draw  a  pure  dielectric  fiber  and  after  drawing  sputter  the  fiber  with  silver  from  one  side,  creating 
split  ring  resonators  (Figure  2).  We  measured  that  these  fibres  have  a  magnetic  resonance  in  the  Thz, 
resulting  in  negative  relative  permeability  [2,3]. 


Figure  2:  Fabrication  of  magnetic  fiber  metamaterial  through  sputtering.  The  silver  coating  is  on  three  sides  of  the 
rectangular  fiber  only. 

The  second  method  was  to  create  a  non-circular  indium  shape  within  a  polymer  preform  and  draw 
this  preform  into  a  fiber.  Figure  3  shows  a  resulting  fiber,  which  also  demonstrated  a  magnetic 
resonance  in  the  THz  (right  panel  of  Fig.  3  shows  the  enhanced  and  phase  shifted  magnetic  field  at 
resonance  [3,4].  At  the  time  of  writing  we  are  exploring  advanced  drawing  techniques  suing  multiple 
polymers  to  improve  the  shape  of  the  magnetic  resonators,  with  extremely  promising  results  that 
indicate  high  quality  magnetic  resonators  should  be  feasible  also  for  the  infrared  spectrum. 
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Figure  3:  Direct  drawn  magnetic  resonator  (left);  transmission  of  an  array  of  the  fiber  and  magnetic  field  at  resonance 
(right). 

We  have  experimentally  demonstrated  that  the  technique  can  draw  materials  with  feature  sizes  well 
below  the  micron  [4],  paving  the  way  to  mid-infrared  and  visible  metamaterials.  We  have 
numerically  demonstrated  that  this  technique  will  enable  the  fabrication  of  intrinsically  invisible 
fibres  [5,6],  by  creating  composite  materials  that  have  the  same  refractive  index  as  air  (Figure  4). 


Figure  4:  Electric  field  propagating  around  (a)  a  silver  wire,  (b)  a  glass  fiber  (c)  a  glass/silver  nano-composite,  all  of  same 
size,  at  a  wavelength  of  633nm  (red  light).  The  wave  fronts  are  unaffected  by  the  composite  fiber,  demonstrating  the 
fiber  is  effectively  invisible  (scattering  cross  section  reduction  of  more  than  95%). 

Future  Work:  There  are  two  main  directions  for  future  work.  The  first  is  to  make  larger  more 
complex  structures  on  the  current  scale  and  use  these  to  realize  practical  devices.  The  second  is  to 
make  smaller  scale  structures  for  use  in  the  IR  and  visible  spectral  regions. 

Conclusion:  The  project  has  been  extremely  successful,  demonstrating  the  feasibility  of  this  novel 
technique  for  the  practical  fabrication  of  metamaterials.  Having  demonstrated  the  feasibility  for 
simple  structures  at  THz  wavelengths,  it  opens  up  the  possibility  of  practical  metamaterial  devices  in 
the  THz,  and  at  shorter  wavelengths. 
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Drawn  metamaterials  with  plasmonic  response  at  terahertz  frequencies 
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Electromagnetic  metamaterials  attract  much  attention  since  they  can  be  engineered  to  exhibit  optical 
properties  not  found  in  nature.  Their  fabrication,  however,  is  challenging,  especially  in  volume.  We 
introduce  drawing  as  a  means  of  fabricating  metamaterials,  thus  demonstrating  a  terahertz 
metamaterial.  We  codraw  polymethyl-methacrylate  and  indium,  producing  several  meters  of 
metamaterial  with  wire  diameters  down  to  —10  gm,  and  lattice  constants  of  —100  jam.  We 
experimentally  characterize  the  transmission  properties  of  different  samples,  observing  high-pass 
filtering  between  0.3-0.4  THz,  in  good  agreement  with  simulations.  ©  2010  American  Institute  of 
Physics,  [doi:  10. 1063/1 .3428576] 


Metamaterials  hold  great  promise  since  they  can  be  en¬ 
gineered  to  exhibit  electromagnetic  properties  which  are  not 
found  in  nature,  by  tailoring  their  subwavelength  structure. 
Recent  demonstrations  of  exotic  structures  using  metamate¬ 
rials  include  electromagnetic  cloaks,1  negative  refractive  in¬ 
dex  media,2  and  sub-diffraction  limited  lenses.3  Such  de¬ 
vices,  however,  are  challenging  to  fabricate,  especially  in 
large  volumes.  In  this  paper  we  report  an  example  of  a  drawn 
metamaterial  which  has  potential  for  inexpensive  volume 
production.  We  produce  more  than  20  m  of  metamaterial  in 
the  form  of  a  terahertz  (THz)  high-pass  filter,  with  potential 
applications  in  THz  imaging4  and  thermal  detection.5 

The  electric  response  of  our  metamaterial  is  tailored  via 
a  two-dimensional  array  of  thin  metallic  wires  with  sub¬ 
wavelength  diameter  and  pitch.  Most  bulk  metals  have 
plasma  frequencies  in  the  visible  and  ultraviolet.  However,  a 
two-dimensional  metallic  wire  array  of  diameter  d  and  sepa¬ 
ration  a  in  a  dielectric  with  permittivity  ed  and  refractive 
index  nd,  irradiated  by  an  electromagnetic  field  with  the 
electric  component  parallel  to  the  wires,  can  be  seen  as  a 
metamaterial  with  lowered  effective  plasma  frequency  fp 
such  that 


f2  _  (C(/«rf)2 

p  27iAcell[ln  a2/d(2a  -  d)]  ’ 


(1) 


where  c0  is  the  speed  of  light  in  vacuum  and  Acell  is  the 
average  area  per  wire  of  the  metamaterial.6  The  permittivity 
as  a  function  of  frequency  /  is  then  given  by 


*eff  if)  = 


f  +  2iyf] 


(2) 


where  y  is  a  damping  factor.6  This  equation  describes  the 
effective  medium’s  transition  from  metal-like  behavior  (low 
transmission  due  to  negative  permittivity),  to  dielectric-like 
behavior  (higher  transmission  due  to  positive  permittivity). 
Arrays  of  thin  metallic  wires  are  thus  metamaterial  high-pass 
filters  in  the  vicinity  of  the  plasma  frequency  fp. 

The  fabrication  of  metamaterials  operating  at  frequen¬ 
cies  higher  than  the  microwave  remains  challenging,  often 
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relying  on  expensive  lithographic  techniques  which  can  only 
be  used  up  to  centimeter  scales.7  We  present  a  fabrication 
approach  inspired  by  fiber  drawing8  used  for  the  large  scale 
and  accurate  fabrication  of  photonic  crystal  fibers  (PCFs). 
Figure  1  (a)  shows  a  schematic  of  how  our  metamaterials  are 
drawn;  a  macroscopically  sized  metal-dielectric  preform  is 
heated  and  reduced  in  size  by  several  orders  of  magnitude. 
However,  there  are  additional  fabrication  challenges  when 
compared  to  conventional  PCFs;  the  scales  required  are  typi¬ 
cally  smaller  (^200  nm  for  optical  metamaterials),  and  be¬ 
ing  metal-dielectric  structures,  material  properties  (such  as 
viscosity  and  melting  temperature)  are  not  well  matched.  In 
spite  of  this,  many  meters  of  microstructured  copper- silica 
fibers  were  recently  drawn  for  the  purpose  of  demonstrating 
guided  modes  via  metallic  reflection.9  Similar  structures 
have  been  achieved  by  pumping  liquid  metals  into  already 
drawn  holey  PCFs  with  submicrometer  holes,10  however  di¬ 
rect  pumping  is  limited  by  pressure  to  fabricating  samples  a 
few  meters  in  length  at  most,  while  direct  drawing  can  in 
principle  be  scaled  to  kilometer  lengths. 

Here  we  present  the  fabrication  of  drawn  metamaterials, 
comprised  of  microstructured  indium-filled  polymethyl¬ 
methacrylate  (PMMA)  fiber,  exhibiting  plasmonic  response 
in  the  THz.  Indium  wires  below  —10  jam  diameter  and 
—  100  jam  separation  are  achieved,  though  submicron  fea¬ 
tures  are  possible.10  We  show  that  the  plasmonic  response 


FIG.  1.  (Color  online)  (a)  Schematic  of  a  metal-dielectric  preform,  drawn 
into  a  metamaterial  via  heating,  (b)  SEM  micrograph  of  a  fabricated 
590  pm  indium-filled  PMMA  fiber  cross-section. 


(b) 
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FIG.  2.  Schematic  of  the  experimental  setup  for  measuring  metamaterial 
fiber  transmittance. 

can  be  controlled  via  the  drawing  speed,  and  compare  results 
with  numerical  simulations. 

Our  fabrication  procedure  is  based  on  adapting  the 
Taylor- wire  process11  to  indium  wires  in  PMMA  tubes,  in 
combination  with  the  stack-and-draw  procedure  used  for  fab¬ 
ricating  PCFs.  PMMA  was  chosen  because  it  is  an  easily 
machined  and  drawn  dielectric;  indium  because  it  has  the 
appropriate  melting  temperature  (156.6  °C)  for  codrawing 
with  PMMA.  Together  their  properties  are  suitable  for  the 
THz  region.  The  metamaterial  fabrication  procedure  was  as 
follows.  We  included  an  indium  wire  (1  mm  diameter, 
99.99%  purity)  in  a  PMMA  tube  with  1.2  mm  inside  diam¬ 
eter  (ID)  and  12  mm  outside  diameter  (OD),  sealed  at  one 
end.  The  tube  was  evacuated  and  drawn  down  at  a  tempera¬ 
ture  of  185  °C  to  an  outer  diameter  of  1.2  mm.  The  resulting 
fiber  was  then  cut  into  smaller  pieces,  which  were  stacked  in 
a  PMMA  capillary  (6.0  mm  ID,  7.5  mm  OD).  This  stack  was 
sealed  at  one  end  and  fed  into  a  furnace  while  being  rapidly 
drawn,  reducing  the  diameter  of  the  structure  by  a  factor  of 
10,  and  increasing  the  length  by  a  factor  of  100.  By  varying 
the  draw  speed,  we  fabricated  fibers  with  diameters  of  590, 
690,  and  770  yam.  The  metal  wires’  diameters  and  their 
separation  scaled  down  correspondingly.  Figure  1  (b)  shows  a 
scanning  electron  microscope  (SEM)  image  of  a  590  yam 
indium-filled  metamaterial  fiber  cross  section.  The  wires 
have  a  lattice  constant  a  —  100  yam  and  diameter  d  —  8  yam. 

We  characterized  the  optical  properties  of  our  samples 
via  THz  time  domain  spectroscopy.12  A  schematic  of  the  ex¬ 
perimental  setup  is  shown  in  Fig.  2.  The  THz  source  was  an 
InGaAs  surface  emitter  pumped  by  80  fs,  800  nm  pulses 
from  a  mode-locked  Ti: Sapphire  laser,  while  the  detector  was 
a  photoconductive  antenna  which  was  gated  by  a  time- 
delayed  pulse  from  the  same  laser.  The  gated  antenna  current 
which  is  proportional  to  the  incident  THz  electric  field  is 
measured  with  a  lock-in  amplifier.  The  entire  THz  transmis¬ 
sion  spectrum  (between  0. 1-1.1  THz)  was  obtained  from  the 
Fourier  transform  of  the  time-delay  scans  between  pump  and 
probe  pulses.  We  cut  the  final  fibers  in  many  3  cm  pieces, 
aligned  them  in  a  one-dimensional  array,  and  measured  the 
transmittance  for  polarized  light  either  parallel  or  perpen¬ 
dicular  to  the  wires.  The  final  sample  sizes  had  areas  of  3 
X  3  cm2,  and  were  attached  to  a  metal  plate  with  a  circular 
aperture  along  the  collimated  THz  beam.  From  repeated 
measurements  we  estimate  a  noise  level  in  transmittance  of 
about  10%  for  frequencies  higher  than  0.2  THz,  and  20%  for 
frequencies  lower  than  0.2  THz. 

Figure  3  (blue  solid  line)  shows  the  measured  transmit¬ 
tance  of  our  590  yam  fiber  for  wires  parallel  to  the  electric 
field  E ,  where  they  form  a  metamaterial.  Very  good  agree- 


FIG.  3.  (Color  online)  Experimentally  measured  metamaterial  transmit¬ 
tance,  normalized  to  transmission  through  aperture  only.  The  electric  field  E 
is  either  parallel  (solid  line)  or  perpendicular  (dashed  line)  to  the  wires.  A 
measurement  for  an  array  of  500  /mm  PMMA  fibers  is  also  shown  (dotted 
line). 

ment  is  found  with  previous  work.13  We  observe  that  trans¬ 
mission  is  extremely  small  at  low  frequencies,  up  to  a  tran¬ 
sition  region  at  —0.3  THz.  At  higher  frequencies  the  sample 
behaves  like  a  dielectric  until  the  wavelength  becomes  too 
short  and  our  structure  can  no  longer  be  seen  as  a  homoge¬ 
neous  medium,  as  discussed  in  detail  below. 

The  dashed  red  line  in  Fig.  3  shows  the  measured  trans¬ 
mittance  when  the  wires  are  perpendicular  to  E.  In  this  case, 
Eqs.  (1)  and  (2)  do  not  apply,  and  the  transmittance  is  close 
to  that  of  an  array  of  500  yam  pure  PMMA  fibers  (black 
dotted  line).  Thus,  our  system  is  both  a  polarizer  and  a  high- 
pass  filter  metamaterial. 

To  study  the  ability  to  engineer  plasmonic  responses  by 
varying  the  draw  speed,  we  measured  the  transmittance  of 
arrays  of  different  fiber  diameters,  with  E  parallel  to  the 
wires.  The  results  are  shown  in  Fig.  4(a).  The  plasmonic 
transition  region  shifts  to  lower  frequencies  as  we  increase 
the  fiber  diameter,  since  the  wires’  separation  a  and  diameter 
d  scale  correspondingly,  leading  to  a  decrease  in  the  plasma 
frequency  of  the  wire  array,  c.f.  Eq.  (1).  We  modeled  the 
transmittance  using  the  finite  element  package  COMSOL, 
mapping  the  micro  structured  geometries  and  wire  diameters 
obtained  from  SEM  images,  and  using  PMMA  refractive 
index14,15  and  indium  permittivity16  [Fig.  4(b)],  observing 
good  agreement  with  the  experiment.  Note  that  Eq.  (1)  pre¬ 
dicts  plasma  frequencies  of  0.59  THz,  0.50  THz,  and  0.45 
THz  for  the  590  yam,  690  yam,  and  770  yam  fibers,  respec¬ 
tively,  assuming  Acell=  ^3a2/2.  Figure  4(a)  shows  that  the 
measured  transmission  peaks  below  such  values.  This  is  due 
to  the  PMMA  jacket  that  surrounds  the  metamaterial  present 
at  the  fiber  core.  To  illustrate  this,  we  simulate  the  transmit¬ 
tance  of  homogeneous  metamaterials  with  plasma  frequen¬ 
cies  fp  predicted  by  Eq.  (1),  possessing  permittivity  de¬ 
scribed  by  Eq.  (2),  surrounded  by  a  PMMA  jacket  [Fig. 
4(c)].  Figures  4(b)  and  4(c)  are  in  excellent  agreement  at  low 
frequencies  (white  background),  where  the  wavelength  is 
much  larger  than  any  feature  of  the  wire  array.  Some  discrep¬ 
ancies  appear  at  higher  frequencies,  when  the  structure’s  fea¬ 
tures  are  no  longer  small  compared  to  the  wavelength 
(shaded  region):  at  these  frequencies  Bragg  scattering  and 
other  resonant  phenomena  due  to  details  of  the  structure  can 
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FIG.  4.  (Color  online)  (a)  Experimental  and  (b)  simulated  transmittance  for 
different  arrays  of  metamaterial  fiber,  (c)  Simulated  transmittance  for 
metamaterials  surrounded  by  a  PMMA  jacket.  The  electric  field  is  parallel  to 
the  wires.  The  shaded  region  illustrates  where  the  medium  cannot  be  seen  as 
homogeneous. 


alter  the  transmission,  and  the  array  can  no  longer  be  con¬ 
sidered  as  a  homogeneous  medium.  The  shaded  region  is  to 
aid  the  eye  and  does  not  represent  a  theoretical  value  at 
which  the  medium  ceases  to  be  homogeneous. 

Note  that  a  similar  metamaterial  was  fabricated  some 
years  ago  using  microstereolithography.  While  such  a 
structure  was  only  of  the  order  of  a  few  millimeters  in  size, 
we  were  able  to  inexpensively  fabricate  many  meters  of  a 
similar  geometry  through  drawing.  Our  technique  potentially 
allows  to  fabricate  kilometers  of  continuous  metamaterial 
with  feature  sizes  down  to  the  nanoscale,  moving  device  op¬ 
erating  frequencies  into  the  mid-infrared  and  visible.  Such 
fiber- shaped  metamaterials  can  be  cut  and  stacked  into  bulk 
materials  to  realize  a  large  number  of  devices,  or  woven  into 
more  complex  structures.  This  could  open  the  way  to  fabri¬ 


cating  cost-effective  plasmonic  devices  and  optical  metama¬ 
terials  on  a  mass  scale,  though  in  this  frequency  range  using 
low  loss  metals  and  dielectrics,  such  as  silver  and  silica, 
would  be  preferable. 

In  conclusion,  we  have  demonstrated  drawn  metamateri¬ 
als  that  act  as  THz  high-pass  filters  and  polarizers.  Transmit¬ 
tance  was  measured  via  THz  time  domain  spectroscopy, 
finding  good  agreement  with  numerical  simulations.  Finally, 
we  have  introduced  drawing  as  a  means  of  mass-producing 
metamaterials.  An  extension  of  this  procedure  would  reduce 
metal  features  down  to  the  nanoscale,  potentially  leading  to 
low-cost  optical  metamaterial  and  plasmonic  devices. 

This  work  was  performed  in  part  at  the  Optofab  node  of 
the  Australian  National  Fabrication  Facility  using  Common¬ 
wealth  and  NSW  State  Government  funding.  This  material 
is  based  on  research  sponsored  by  the  Air  Force  Research 
Laboratory,  under  Agreement  No.  FA23 86-09- 1-4084. 
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Abstract:  We  present  the  design  of  an  invisible  metamaterial  fibre 
operating  at  optical  frequencies,  which  could  be  fabricated  by  adapting 
existing  fibre  drawing  techniques.  The  invisibility  is  realised  by  matching 
the  refractive  index  of  the  metamaterial  fibre  with  the  surroundings.  We 
present  a  general  recipe  for  the  fabrication  of  such  fibres,  and  numerically 
characterise  a  specific  example  using  hexagonally  arranged  silver  nanowires 
in  a  silica  background.  We  find  that  invisibility  is  highly  sensitive  to  details 
of  the  metamaterial  boundary,  a  problem  that  is  likely  to  affect  most 
invisibility  and  cloaking  schemes. 
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1.  Introduction 

Electromagnetic  invisibility,  i.e.  the  passing  of  light  through  a  structure  without  scattering, 
distortion,  or  absorption,  is  currently  of  great  interest,  due  to  its  fascinating  implications.  A 
number  of  structures  that  are  invisible  to  electromagnetic  radiation  have  been  theoretically 
proposed  and  experimentally  demonstrated  in  the  past  few  years;  from  simple  devices  such  as 
antireflection  coatings  [1]  and  absorbing  screens  [2],  to  more  exotic  structures  such 
electromagnetic  cloaks,  which  can  hide  an  object  with  a  plasmonic  medium  [3-5]  or  a 
metamaterial  possessing  especially  engineered  effective  optical  properties  [6-10].  In  this 
paper  we  discuss  metamaterials  that  are  themselves  optically  transparent,  as  opposed  to 
possessing  cloaking  properties.  Recent  proposals  of  transparent  metamaterials  are  based  on 
alternating  layers  of  natural  materials  and  double-negative  metamaterials  [11,12]  and 
transformation  optics  [13]. 

Here  we  propose  a  much  simpler  approach:  optical  transparency  is  achieved  by  designing 
a  metamaterial  with  refractive  index  that  matches  the  surroundings,  consisting  of  an  array  of 
sub- wavelength  metallic  nanocylinders  surrounded  by  a  dielectric.  Such  structures  could 
potentially  be  fabricated  via  metamaterial  fibre  drawing,  either  via  direct  co-drawing,  in 
which  a  macroscopically  sized  metal-dielectric  preform  is  heated  and  reduced  in  size  by 
several  orders  of  magnitude  [14,15],  or  via  pumping  liquid  metal  into  existing  micro-  and 
nano-  structured  holey  photonic  crystal  fibres  (PCFs)  [16].  Holey  fibres  containing  feature 
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Fig.  1.  Schematic  of  a  metal-dielectric  preform,  drawn  into  a  metamaterial  via  heating, 
designed  to  be  invisible  at  a  chosen  optical  wavelength  (e.g.  633nm). 

sizes  down  to  30nm  have  recently  been  reported  [17],  showing  that  drawn  patterned 
nanostructures,  and  thus  drawn  optical  metamaterials,  can  in  principle  be  produced. 

Fibre  drawing  has  already  proven  extremely  successful  for  the  large  scale  and  accurate 
fabrication  of  PCFs  [18].  In  the  context  of  metamaterials,  drawing  has  significant  advantages 
over  conventional  fabrication  techniques,  especially  in  terms  of  cost  (the  macroscopic  preform 
is  generally  an  inexpensive  dielectric  a  few  cm3  in  volume)  and  volume  production  (one 
preform  can  in  principle  be  used  to  make  kilometers  of  fibre).  In  contrast,  the  lithographic 
techniques  used  for  producing  metamaterials  at  operating  frequencies  from  the  terahertz  to  the 
visible  are  generally  quite  expensive,  and  produce  no  more  than  a  few  cm2  of  metamaterial  at 
a  time  [19].  However,  while  drawing  offers  the  possibility  of  creating  nanoscale  structures  in 
volume,  this  comes  with  a  constraint  -  they  have  to  be  uniform  along  one  dimension.  Recently 
it  has  been  shown  that  important  functionality  can  potentially  be  realised  with  such  fibre- 
based  metamaterials,  such  as  low-loss  mid-IR  waveguiding  [20],  sub-wavelength 
waveguiding  [21]  and  hyperlensing  [22]. 

Here  we  propose  another  interesting  structure  that  can  be  fabricated  via  drawing  of  metal 
cylinders  in  a  dielectric  -  the  optically  invisible  fibre.  Figure  1  illustrates  a  schematic  of  the 
potential  fabrication  process,  whereby  a  macroscopic  metal-dielectric  preform  is  drawn  into  a 
fibre  that  is  invisible  around  a  chosen  optical  wavelength  (e.g.  X  =  633nm  for  a  He:Ne  laser) 
with  normal  incidence.  Our  design  method  is  based  on  matching  the  effective  optical 
parameters  for  a  metamaterial  of  2-dimensional  patterned  nanocylinders  of  sub-wavelength 
pitch  to  its  surroundings.  We  characterise  our  invisible  fibre  via  the  total  scattering  cross 
section  oT ,  showing  that  cT  has  a  minimum  when  the  real  part  of  the  metamaterial  effective 
refractive  index  matches  the  surroundings.  Taking  the  example  of  a  fibre  to  be  invisible  at 
633nm  using  silver  nanocylinders  in  a  silica  cylinder  of  radius  R  =  1pm,  we  study  the 
behaviour  of  oT  with  wavelength,  incident  angle,  and  polarisation.  In  this  case  more  than  a 
95%  reduction  in  oT  can  be  achieved  over  ~13nm  bandwidth  when  compared  to  dielectric  and 
conducting  cylinders  of  the  same  size,  as  well  as  a  95%  reduction  in  oT  for  angles  that  are  ± 
10°  to  normal  incidence  at  633nm.  We  find  that  invisibility  is  highly  sensitive  to  changes  in 
the  metamaterial  boundary,  with  best  results  when  the  total  filling  fraction  -  including  the 
metamaterial  boundaries  -  matches  the  bulk  filling  fraction.  Thus,  particular  attention  should 
be  paid  to  the  fibre  preform  design  process.  At  larger  fibre  diameters  our  invisible 
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metamaterial  becomes  a  near-perfect  absorber,  due  its  low  reflectance  and  inherently  lossy 
nature. 

2.  Material  design  -  bulk  properties 

2.1  Principle 

An  object  in  vacuum  is  invisible  if  it  neither  reflects  nor  absorbs  light,  and  if  it  induces  no 
changes  in  phase  with  respect  to  its  surroundings.  Our  objective  is  to  realise  an  optical 
metamaterial  that  closely  meets  such  requirements. 

The  optical  properties  of  a  homogeneous  material  are  commonly  characterised  via  its 
complex  impedance  z  and  refractive  index  n ,  related  to  the  relative  electric  permittivity  s  =  n/z 
and  relative  magnetic  permeability  p  =  nz.  The  impedance  mismatch  at  an  interface  causes 
reflections,  while  the  refractive  index  characterises  the  phase  velocity  in  the  medium.  For 
example,  for  a  homogeneous  slab  of  length  L,  the  transmission  and  reflection  coefficients  t 
and  r  for  a  plane  wave  of  angular  frequency  co  at  normal  incidence  are  given  by  [23] 


_ l _ 

ncoL  if  1  ^  .  ncoL  ’ 

cos - z  +  —  sin - 

c  2\  z)  c 


(1) 


(2) 


where  c  is  the  speed  of  light  in  vacuum. 

It  goes  without  saying  that  a  hypothetical  slab  with  n=  1  and  z—  1  is  invisible  in  vacuum, 
since  r  -  0  and  t  =  exp(icolVc).  Natural  materials  possess  an  impedance  and  refractive  index 
that  are  not  unity  at  optical  wavelengths,  and  are  thus  visible;  typically  bulk  metals  possess  a 
permittivity  sm<  0  and  dielectrics  possess  sd>  1 ,  with  p  =  1  for  both. 

However,  it  is  well  known  that  a  composite  medium  composed  of  a  2-dimensional  array  of 
thin  metal  wires  surrounded  by  a  dielectric  in  the  electrostatic  limit  (feature  sizes  much 
smaller  than  the  wavelength)  irradiated  by  TE  polarised  light  (electric  field  directed  along  the 
wires)  can  be  treated  as  a  homogeneous  medium  with  effective  permittivity  £eff  given  by 

Seff=fsm+{\~f)sd,  (3) 


where/is  the  area  filling  fraction  of  the  metal  (see,  for  example,  Ref  [24].).  Thus,  by  knowing 
sm  and  sd  at  a  given  wavelength,  it  is  possible  to  compute  the  filling  fraction  at  which  the 
composite  medium  is  invisible  by  solving  for  seff=  1  in  Eq.  (3),  yielding 


/=  — —■  (4) 

£  —£. 

m  d 

Assuming  for  now  that  the  composite  medium  is  effectively  non-magnetic  (peff  =  1),  for  / 
satisfying  Eq.  (4)  the  effective  index  and  impedance  take  the  values  neff  =  1  and  zeff  =  1 
respectively,  as  required. 

By  applying  the  above  procedure  at  a  wavelength  A  =  633nm  using  silver 
( sAg  =  -11. 7  +  0.49/  from  interpolated  values  of  Ref  [25].)  and  silica  (ssmca  =  2.123  from  the 
Sellmeier  equation  [26]),  one  obtains  that  £eff  =  1.00  +  0.03/  when /=  0.0567.  For  a  hexagonal 

lattice  of  cylinders  with  radius  r  and  centre-to-centre  pitch  d ,  i.e.  /  =  nr1  /  (V3 d2  /  2) , 
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0 


Fig.  2.  (a)  Real  and  (b)  imaginary  parts  of  the  permittivity  of  silver  at  optical  wavelengths:  the 
blue  solid  line  represents  the  polynomial  fit  used  in  the  numerical  model,  the  red  markers 
indicate  experimental  values. 


one  finds  that  for  a  cylinder  diameter  of  2 r  =  30nm  (the  smallest  drawn  feature  reported  [17]), 
the  pitch  must  be  d  =  120nm,  in  which  case  XI d  ~4  and  the  electrostatic  approximation  is  not 
appropriate. 

Thus,  while  Eq.  (3)  is  a  good  starting  point,  higher  order  corrections  [27]  are  necessary  to 
achieve  higher  precision  in  calculating  the  filling  fraction  required  for  invisibility. 

2.2  Higher  order  corrections 

Here  we  obtain  corrections  to  the  first-order  theory  by  directly  calculating  r  and  t  and 
inverting  Eqs.  (1)  and  (2), 

ncoL  1  rl  2  2i 

cos - =  — [1-r  +r],  (5) 

c  2 1 


l(l  +  r)2-t 

Ala -rf-t2  ’ 


for  different  slabs  of  lengths  L  =  d,  2d,  3d.  Following  the  procedure  outlined  in  [23],  using 
Eq.  (5)  and  (6)  we  verify  that  z  is  independent  of  L,  and  unambiguously  retrieve  n  (and  thus  s 
and  n)  at  a  given  wavelength. 

We  use  a  multipole  expansion  in  combination  with  the  transfer  matrix  method  to  calculate 
r  and  t  for  uniform  slabs  of  hexagonally  arranged  silver  cylinders  in  a  silica  background  under 
TE  polarisation.  Fields  around  each  cylinder  are  expanded  in  a  multipole  basis,  where  analytic 
expressions  are  derived  for  the  cylindrical  boundary  conditions  [28].  Lattice  sums  are  used  to 
relate  the  field  incoming  onto  each  cylinder  to  the  field  outgoing  from  other  cylinders  within  a 
singer  layer  of  the  lattice  [29].  The  scattering  and  reflection  coefficients  of  a  multi-layer  stack 
of  gratings  are  computed  from  the  coefficients  of  a  single  grating  layer  via  a  recurrence 
relation.  The  permittivity  of  silica  is  included  in  the  calculations  via  a  wavelength-dependent 
Sellmeier  expansion  [26],  whereas  for  silver  we  use  a  ninth-order  polynomial  fit  to  measured 
real  and  imaginary  parts  of  complex  permittivity  [25]  [Fig.  2(a)  and  2(b)]. 

We  use  this  method  to  calculate  r,  t,  and  extract  n ,  z,  as  a  function  of/ for  two  fixed  unit 
cell  lengths,  d  =  lOnm  and  d  =  lOOnm.  The  results  are  shown  in  Fig.  3(a)  and  3(b):  whereas 
for  smaller  lattices  ( d  =  lOnm,  i.e.  XI d  ~60),  calculations  are  in  excellent  agreement  with  the 
quasistatic  theory  of  Eq.  (3),  for  larger  lattices  ( d  =  lOOnm,  i.e.  X/d  ~6)  a  significant  correction 
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Fig.  3.  (a)  Real  and  (b)  imaginary  parts  of  effective  refractive  index  and  impedance  as  a 
function  of  filling  fraction /for  hexagonally  arranged  silver  cylinders  in  a  silica  background  at 
a  wavelength  of  X  =  633nm,  retrieved  for  d  =  lOnm  and  d  =  lOOnm,  and  calculated  through 
Eq.  (3).  Note  that  effective  medium  theory  [Eq.  (3)]  applies  only  for  smaller  unit  cells,  (c)  Real 
parts  of  retrieved  effective  optical  parameters  as  a  function  of  wavelength  for  d  =  lOnm,  /  = 

0.0567  and  (d)  for  d  =  lOOnm,  /=  0.0659.  In  both  cases  Re(n)  =  1  at  633nm  by  design,  but  for 
the  larger  unit  cell  size  a  magnetic  response  is  present. 

to  the  first  order  theory  is  required.  Analytical  work  on  the  homogenisation  of  arrays  of 
metallic  cylinders  [27]  has  shown  that  lowest  order  dynamic  corrections  to  the  static 
homogenisation  theory  lead  to  an  additional  dipolar  magnetic  response  as  well  as  a 
quadrupolar  electric  response,  both  of  same  order  of  magnitude.  In  particular,  the  combination 
of  these  two  effects  explains  that  we  have  z  ^  1  In. 

From  Fig.  3(b)  we  conclude  that  a  filling  fraction  /  =  0.0659  yields  an  effective  index  of 
1.00  +  0.01/  at  X  =  633nm  for  a  lattice  constant  of  d  -  lOOnm,  corresponding  to  a  cylinder 
diameter  of  2 r  =  27nm,  which  can  potentially  be  achieved  using  fibre  drawing  techniques. 

Figures  3(c)  and  3(d)  shows  the  real  parts  of  the  optical  parameters  extracted  as  a  function 
of  wavelength  for  d  =  lOnm,  /=  0.0567  and  for  d  =  lOOnm,  /=  0.0659.  In  the  former  case 
[Fig.  3(c)],  all  optical  parameters  are  unity  at  633nm  by  design.  In  the  latter  case  [Fig.  3(d)], 
the  presence  of  a  magnetic  response  makes  it  impossible  to  simultaneously  match  the 
refractive  index  and  the  impedance  to  the  surroundings.  As  we  will  see  next,  this  does  not 
compromise  device  performance  at  the  desired  wavelength.  We  now  characterise  its  visibility 
in  detail  in  a  fibre  geometry. 
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Fig.  4.  (a)  Schematic  of  incident  wave  parameters  with  respect  to  the  axis  of  a  dielectric 
cylinder  with  radius  R,  filled  with  hexagonally  arranged  metal  nanocylinders  of  radius  r  and 
centre-to-centre  pitch  d.  (p  indicates  azimuthal  dependence,  6  indicates  the  incident  angle  of  the 
wave  vector  k  with  respect  to  the  cylinder  axis  and  S  indicates  the  angle  of  polarisation,  (b) 
Normalised  scattering  cross  section  for  different  cylinders  with  R  =  1pm.  The  designed 
metamaterial  cylinder  has  a  minimum  cross  section  at  the  chosen  wavelength  of  633nm.  Inset: 
Normalised  scattering  cross  section  (full  calculation)  for  the  metamaterial  cylinder  as  a 
function  of  azimuthal  angle. 


3.  Numerical  simulations  and  analysis 

We  use  the  multipole  method  [30,31]  to  evaluate  scattering  of  waves  [32]  at  optical 
wavelengths  for  two-dimensional  cylinders  of  radius  R  made  out  of  the  metamaterial 
characterised  in  Section  2.2.  This  method  allows  to  characterise  visibility  via  semi-analytical 
calculations  of  the  scattering  cross  section  (SCS)  for  any  linear  polarisation  state,  azimuthal 
angle,  and  incident  angle  with  respect  to  the  cylinder  axis  [represented  by  S ,  (p  and  6 
respectively,  as  shown  in  the  schematic  of  Fig.  4(a)].  The  scattering  cross  section  a  is  defined 
as  the  total  scattered  power  per  unit  length  at  large  distances,  normalised  to  the  incident  power 
density.  In  polar  coordinates  (p,a)  this  is  written  as 

\Es(p,a)\2 

a  (a)  =  lim  2 np - — ,  (7) 

|E°(p,a)| 

where  E°  and  Es  are  the  incident  and  scattered  electric  field,  respectively.  One  then  obtains 
the  total  SCS  oT  by  integrating  over  all  angles, 


2  n 


LK 

|  cr(a)da . 


(8) 


Another  parameter  for  characterising  the  scattering  properties  of  objects  is  the  equivalent 
blockage  width  Weq  a  complex-valued  parameter  defined  as  the  width  of  an  ideal  shadow 
which  produces  the  same  forward-scattered  field  as  the  cylinder  being  observed  [33].  It  can  be 
shown  that  oT  =  2 Re(Weq).  It  follows,  for  example,  that  at  a  given  wavelength  a  perfectly 
conducting  cylinder  of  radius  R  possesses  (jj/2R  =  2  in  the  limit  of  R^> oo,  whereas  a  perfectly 
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Fig.  5.  Electric  field  (out  of  plane)  at  X  =  633nm  incident  on  a  R  =  1pm  (a)  silver,  (b)  silica  and 
(c)  metamaterial  cylinder,  (d)  Magnitude  of  scattered  electric  field  for  1pm  metamaterial 
cylinders  at  A  =  593nm,  where  only  the  impedance  is  unity,  and  (e)  at  2  =  633nm,  where  only 
the  refractive  index  is  unity. 


absorbing  cylinder  possesses  gtI2R  =  1  in  the  same  limit.  We  refer  to  the  vast  literature  on  the 
scattering  cross  section  for  further  details  [34] . 

We  calculate  oT  at  optical  wavelengths  for  a  fibre  composed  of  the  metamaterial  presented 
in  Section  2.2,  formed  by  361  hexagonally  arranged  silver  cylinders  with  radius  r  =  13.5nm 
and  pitch  d  =  lOOnm  within  a  silica  cylinder  of  radius  R  =  1pm  in  vacuum.  We  consider 
normal  incidence,  with  the  electric  field  parallel  to  the  cylinders  (TE  polarisation,  6  =  90°,  S  = 
90°).  The  results  are  shown  in  Fig.  4(b):  the  metamaterial  fibre  has  gtI2R  =  -12.6dB  at  2  = 
633nm,  corresponding  to  more  than  a  97%  reduction  in  shadow  when  compared  to  equivalent 
silver  or  silica  cylinders  of  the  same  size.  Furthermore,  the  metamaterial  fibre  has  gtI2R  <0.1 
over  a  13nm  bandwidth  (95%  reduction  in  shadow).  Additionally,  an  effective  reduction  of  gt 
occurs  over  -  lOOnm  with  respect  to  either  of  its  constituent  materials.  The  inset  of  Fig.  4(b) 
presents  the  normalised  total  SCS  at  normal  incidence  for  2  =  633nm,  as  a  function  of  the 
azimuthal  angle  cp,  demonstrating  -12  +  0.6  dB  reduction  of  total  SCS  at  all  angles.  The 
corresponding  contour  plot  of  the  electric  field  at  2  =  633nm  for  these  three  cylinders  is  shown 
in  Fig.  5(a)-5(c);  the  designed  fibre  indeed  appears  invisible,  exhibiting  low  scattering  and 
loss. 

We  also  calculate  gt  for  the  same  cylinder  possessing  an  effective  s  and  //  retrieved  in 
Section  2.2  [Fig.  3(d)],  finding  excellent  agreement  with  the  full  multipole  calculation 
[Fig.  4(b)].  The  minimum  SCS  occurs  at  a  wavelength  when  Re(n)  =  1  (minimum  phase 
change),  as  opposed  to  when  Re(z)  =  1  (minimum  reflection).  Thus,  the  interference  due  to  a 
phase  change  through  the  fibre,  due  to  a  small  perturbation  in  n ,  has  a  more  significant  effect 
on  gt  than  a  reflection  due  to  a  small  perturbation  in  z.  This  is  further  illustrated  in  Fig.  5(d) 
and  5(e),  which  compares  the  scattered  electric  field  magnitude  IESI  for  the  metamaterial  fibre 
at  2  =  593nm  (Re(z)  =  1.00,  Re(n)  =  1.08)  and  2  =  633nm  (Re(z)  =  1.07,  Re(n)  =  1.00). 
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Fig.  6.  (a)  Colour  density  plot  of  normalised  SCS  as  a  function  of  incident  angle  0  for  TE 
polarisation  (8  =  90°),  and  (b)  as  a  function  of  polarisation  angle  S  for  normal  incidence  (0  = 
90°). 


4.  Incident  angle,  polarisation,  and  size  dependence 

We  now  examine  the  dependence  of  gt  on  polarisation  and  incident  angle  with  respect  to  the 
cylinder  axis  (S  and  6  respectively).  Thus  far  we  have  considered  the  case  3  =  90°,  9  =  90°, 
corresponding  to  TE  polarisation  at  normal  incidence.  Figure  6(a)  shows  a  density  plot  of  total 
scattering  cross  section  as  a  function  of  incident  angle,  while  maintaining  the  polarisation  (i.e. 
we  fix  3  =  90°  and  vary  0).  Note  that  gt  <  -lOdB  for  9  =  90  +  10°  with  respect  to  the  normal 
at  X  -  633nm.  As  we  increase  the  incident  angle,  the  regions  of  minimum  scattering  move  to 
longer  wavelengths,  and  for  6  =  90  +  50°  there  exists  a  wavelength  range  where  gtI2R  <  1 
over  a  bandwidth  of  ~10nm. 

Figure  6(b)  shows  a  density  plot  of  total  scattering  cross  section  as  a  function  of 
polarisation,  at  constant  incident  angle  (i.e.  we  fix  0  =  90°  and  vary  8).  As  expected,  the 
regions  of  smallest  scattering  correspond  to  TE  polarisation.  As  the  electric  field  along  the 
fibre  decreases,  so  does  the  gt.  For  purely  TM  polarisation  (magnetic  field  directed  along  the 
wires),  the  cylinder  behaves  like  a  dielectric. 

In  Section  2  we  found  that  the  metamaterial  used  to  design  our  invisible  fibre  possesses  a 
refractive  index  of  n  =  1.00  +  0.01/  at  2  =  633nm,  corresponding  to  a  loss  of  0.43dB/pm; 
achieving  a  low-loss  transparent  metamaterial  is  thus  only  truly  possible  at  smaller  fibre  radii. 
We  now  examine  the  influence  of  metamaterial  fibre  dimensions  on  the  total  scattering  cross 
section. 

Figure  7(a)  shows  the  normalised  total  SCS  as  a  function  of  R,  calculated  analytically  at  X 
=  633nm  using  retrieved  optical  parameters  for  the  designed  metamaterial  fibre.  As  we 
increase  cylinder  radius  R ,  gt/2R  reaches  unity,  thus  behaving  like  a  near-perfect  absorber.  In 
contrast,  the  normalised  SCS  for  silver  and  silica  cylinders  tends  to  2  at  large  radii,  as 
expected. 

Figure  7(b)  (crosses)  shows  the  full  calculation  of  the  normalised  total  SCS  as  a  function 
of  R ,  using  the  multipole  method.  The  fibres  are  silica  cylinders  of  radius  R  centred  in  the 
origin,  containing  silver  nanocylinders  with  r  =  13.5nm  arranged  in  a  hexagonal  array  of  pitch 
d  =  lOOnm,  such  that  nanocylinder  distance  from  the  origin  is  p  <  R-r.  Since  the  number  of 
cylinders  N  increases  stepwise  with  R ,  the  total  filling  fraction  fT  =  Nr  /R2  when  increasing  R 
is  not  constant,  and  goes  through  discontinuities.  When  fT  differs  from  the  ideal  value  of  /, 
regions  at  the  boundary  exist  for  which  the  local  effective  index  near  the  boundary  of  the 
metamaterial  fibre  is  not  unity,  causing  increased  scattering.  Minimal  gt  is  thus  achieved  at 
values  of  R  for  which  fT  =  0.659.  In  particular,  we  observe  the  lowest  normalised  SCS  of 
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Fig.  7.  (a)  Analytically  computed  oT  at  X  =  633nm  as  a  function  of  cylinder  radius,  for  silver, 
silica,  and  our  metamaterial  cylinder  using  retrieved  parameters,  (b)  Comparison  with  full 
multipole  calculation  for  R  <  3  pm.  Inset:  electric  field  (out  of  plane)  for  R  =  470nm  and  R  = 
480nm. 


-24dB  for  a  metamaterial  cylinder  of  R  =  230nm.  Despite  possessing  significantly  lower  SCS 
compared  to  silver  and  silica  cylinders  of  the  same  size,  oT  is  strongly  sensitive  to  total  filling 
fraction,  since  scattering  occurs  when  the  local  filling  fraction  at  the  metamaterial  boundary 
significantly  deviates  from  that  required  by  design.  For  example,  the  inset  of  Fig.  7  shows  a 
contour  plot  of  the  electric  field  for  R  =  470nm  (containing  73  nanocylinders,  thus 
fT  =  0.602  and  aT!2R  =  -6.2dB)  and  R  =  480nm  (containing  85  nanocylinders,  thus  fT  =  0.672 
and  a  much  reduced  normalised  total  SCS,  gtI2R  =  -19.4dB). 

Note  that  the  metamaterial  fibre  with  R  =  3pm  already  contains  3241  silver  nanocylinders. 
Due  to  memory  constraints,  we  could  not  perform  a  full  calculation  of  the  metamaterial  total 
SCS  for  R  >  3  pm. 

6.  Conclusions  and  future  outlook 

In  conclusion,  we  have  designed  and  characterised  a  metamaterial  fibre  which  exhibits  a 
strong  reduction  in  scattering  cross  section  compared  to  metal  and  dielectric  cylinders  of 
equivalent  size.  The  design  procedure  can  be  appropriately  adapted  to  create  invisible  fibres  at 
any  optical  wavelength  by  changing  the  metal/dielectric  filling  fraction,  or  create  fibres  with 
other  arbitrary  permittivity  between  that  of  the  background  and  of  the  metal.  Our  device 
operates  optimally  at  the  designed  wavelength  under  TE  polarisation  and  normal  incidence, 
yet  excellent  performance  is  maintained  for  ~10nm  of  bandwidth,  or  within  a  ±  10°  variation 
in  incident  angle  with  respect  to  normal  incidence,  where  a  reduction  in  shadow  of  more  than 
95%  is  maintained.  Additionally,  visibility  is  significantly  reduced  over  lOOnm  bandwidth 
with  respect  to  its  constituent  materials.  Optical  parameters  obtained  from  extraction 
procedures  give  a  good  estimate  of  the  total  scattering  cross  section,  unless  imperfect 
termination  of  the  nanostructure  at  the  edges  causes  additional  scattering.  As  with  other 
invisible  metamaterials,  our  device  is  intrinsically  limited  by  loss,  being  invisible  at  small 
radii,  yet  behaving  like  a  near-perfect  absorber  for  larger  fiber  radii.  The  feasibility  of  co¬ 
drawing  metal-dielectric  micro  structured  fibres  has  already  been  demonstrated  [14,15],  and  so 
has  the  feasibility  of  drawing  air-dielectric  nanoscopic  features  [17];  future  work  will  focus  on 
combining  these  techniques  to  produce  nanoscale  metallic  features.  While  being  harder  to 
fabricate,  the  concept  we  presented  can  be  extended  to  produce  other  3 -dimensional  optically 
invisible  structures,  such  as  spheres,  which  would  have  a  wider  tolerance  on  incident  angle 
and  polarisation. 
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Abstract:  We  present  a  novel  method  for  producing  metamaterials  with  a 
terahertz  magnetic  response  via  fiber  drawing,  which  can  be  inexpensively 
scaled  up  to  mass  production.  We  draw  a  centimeter  preform  to  fiber,  spool 
it,  and  partially  sputter  it  with  metal  to  produce  extended  slotted  resonators. 
We  characterize  metamaterial  fiber  arrays  with  different  orientations  via 
terahertz  time  domain  spectroscopy,  observing  distinct  magnetic  resonances 
between  0.3  and  0.4  THz,  in  excellent  agreement  with  simulations. 
Numerical  parameters  retrieval  techniques  confirm  that  such  metamaterials 
possess  negative  magnetic  permeability.  Combined  with  fiber-based 
negative  permittivity  materials,  this  will  enable  the  development  of  the  first 
woven  negative  index  materials,  as  well  as  the  fabrication  of  magnetic 
surface  plasmon  waveguides  and  sub  wavelength  waveguides. 
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1.  Introduction 

Metamaterials  are  artificial  materials,  composed  of  sub-wavelength  “atoms”  that  can 
collectively  exhibit  effective  electromagnetic  responses  not  found  in  nature.  A  large  number 
of  extraordinary  metamaterial  devices  have  been  recently  reported,  e.g.  superlenses  [1], 
electromagnetic  cloaks  [2]  and  negative  refractive  index  materials  [3].  Producing  terahertz 
(THz)  and  optical  metamaterials  usually  requires  expensive  and  labor-intensive  micro-  and 
nano-  fabrication  processes,  such  as  electron-beam  lithography,  focused-ion  beam  milling, 
nanoimprint  lithography  and  direct  laser  writing,  which  limit  metamaterial  sample  sizes  to  a 
few  centimeters  [4].  Here  we  present  a  method  for  the  production  of  a  fiber  metamaterial 
possessing  negative  magnetic  permeability,  p,  in  the  THz.  Fiber  drawing  has  emerged  as 
means  of  realizing  inexpensive  metamaterials,  where  a  macroscopic  (centimeter- sized) 
preform  is  heated  until  softened,  and  drawn  into  a  fiber  that  is  reduced  in  area  by  several 
orders  of  magnitude  and  greatly  increased  in  length.  Although  continuous  metal  wires  down 
to  the  micrometer  [5-9]  and  nanometer  [10]  scale  have  been  reported  using  such  techniques, 
giving  the  ability  to  tailor  the  effective  electric  permittivity,  8,  this  approach  has  previously 
never  been  used  for  creating  magnetically  responsive  metamaterials. 

2.  Fabrication  and  sample  preparation 

A  common  approach  to  produce  a  negative  p  material  is  to  use  sub -wavelength  metallic  split 
ring  resonators  (SRRs)  composed  of  a  planar  open  conductive  split  loop  with  an  inductance  L 
and  a  capacitance  C  that  opposes  incident  magnetic  fields  at  frequencies  near  the  LC 
resonance.  The  geometry  considered  here  however  is  the  slotted  resonator.  This  has  a 
geometry  similar  to  that  of  SRRs,  but  instead  of  being  inscribed  in  a  plane,  has  a  finite  depth 
and  is  invariant  along  the  longitudinal  axis.  Slotted  resonators  can  also  exhibit  a  negative 
magnetic  response  [11],  however  their  depth  demands  new  fabrication  strategies.  Recently, 
slotted  resonators  in  the  near-  and  mid-infrared  were  produced  via  a  combination  of  polymer 
direct-laser  writing  and  chemical  vapor  deposition  [12].  Our  approach  is  to  fabricate  slotted 
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Fig.  1 .  Schematic  of  the  fabrication  procedure,  i)  A  dielectric  square  preform  is 
drawn,  ii)  sputtered  with  silver  on  three  sides  and  iii)  spooled  into  an  array,  (b) 
Optical  microscope  image  of  the  metamaterial  array,  (c)  SEM  micrograph  of  i) 
the  cross  section  of  the  metamaterial  fiber  and  ii)  of  its  250nm  silver  coating. 


resonators  in  the  fiber  form,  utilizing  the  longitudinal  uniformity  of  the  fiber  drawing  process 
to  create  kilometers  of  the  required  structure. 

Figure  1(a)  shows  a  schematic  of  our  fabrication  procedure.  A  preform  with  a  square 
cross  section  with  2  cm  side  was  milled  from  a  rod  of  extruded  cyclo-olefin  Zeonex  [13],  a 
drawable  polymer  with  low  loss  in  the  THz  region.  A  square  cross-section  was  chosen  to 
facilitate  the  spooling  of  the  resulting  fiber  with  a  controlled  orientation.  The  preform  was 
drawn  to  125+5  pm  side,  and  spooled  onto  a  cylinder  which  was  subsequently  loaded  in  a  DC 
magnetron  sputtering  deposition  system  to  coat  a  200-250  nm  silver  thin  film  on  three  sides. 

To  ensure  the  strongest  magnetic  resonance  possible  with  this  geometry  [14],  the  silver 
coating  thickness  was  selected  to  be  larger  than  the  skin  depth  8  at  0.1-1  THz.  A  good 
approximation  of  the  skin  depth  for  this  range  of  frequencies  can  be  obtained  using  the  d.c. 
conductivity  of  silver  [15],  adc  =  6xl07  Q'W1,  through  8  =  l/(/7ip0adc)1/2  with  p0  the 
permeability  of  free  space  and /the  frequency  at  which  the  skin  depth  is  defined.  This  yields  a 
skin  depth  between  50  and  200  nm  over  the  frequency  range  of  interest.  To  improve  adhesion, 
a  chromium  interlayer  (~  10-20  nm)  was  applied  prior  to  the  silver  layer  deposition.  Good 
silver  coating  coverage  of  the  three  exposed  sides  was  ensured  by  translating  and  rotating  the 
spool  inside  the  sputtering  chamber  [Figure  1(a)].  The  resulting  fiber  formed  a  single  slotted 
square  prism  resonator  that  is  invariant  along  the  axis  of  the  fiber.  The  coated  fiber  was  then 
spooled  again  to  produce  a  single  flat  array,  which  was  supported  in  a  frame  [Figure  1(b)]. 
Figure  1(c)  shows  a  scanning  electron  microscope  (SEM)  image  of  the  cross  section  of  our 
final  fiber  with  125  pm  square  edge,  and  a  250  nm  thick  silver  coating. 

By  controlling  the  orientation  of  the  coated  fiber  at  the  start  of  the  second  spooling 
process  it  was  possible  to  produce  arrays  with  different  resonator  orientations,  each  with 
different  electromagnetic  properties,  depending  on  whether  the  resonators  are  symmetric  or 
asymmetric  with  respect  to  the  forwards-  and  backwards -propagating  wave  [16].  Figure  2(c) 
illustrates  the  resonator  array  geometries  studied  here.  In  both  cases,  we  consider  fields  under 
normal  incidence  (i.e.  propagating  in  the  x-direction),  with  either  the  magnetic  field  (TM 
polarization)  or  the  electric  field  (TE  polarization)  directed  along  the  fiber  in  the  z-direction. 
In  the  symmetric  case,  the  resonator  will  respond  isotropically  along  the  x-axis  (i.e.  the 
transmission  and  reflection  will  be  identical  whether  illumination  is  from  the  left  or  the  right), 
whereas  in  the  asymmetric  case  it  will  respond  bianisotropically  -  transmission  and  reflection 
will  depend  on  whether  the  light  comes  from  the  left  or  right.  This  is  due  to  the  lack  of 
inversion  symmetry  and  leads  to  a  magnetoelectric  coupling:  the  magnetic  field  in  z  induces 
an  electric  dipole  in  y,  and  the  electric  field  in  y  also  induces  a  magnetic  dipole  in  z,  in 
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Fig.  2.  (a)  Experimental  and  simulated  transmittance  and  phase  (inset)  under 
TM  polarization.  The  transmission  dip  indicates  a  magnetic  resonance,  (b) 
Experimental  and  simulated  transmittance  under  TE  polarization,  (c)  Schematic 
of  the  fiber  array  orientations,  d)  Relative  phase  of  the  simulated  magnetic  field 
at  the  center  of  the  resonator  with  respect  to  the  incident  field,  e)  Color  plot  of 
the  simulated  magnetic  field  incident  under  TM  polarization  for  the  symmetric 
resonator  at  (i)  0.15  THz,  (ii)  0.34  THz  and  (iii)  0.37  THz. 


particular  leading  to  two  different  characteristic  impedances  for  two  waves  propagating  in 
opposite  directions  [16].  This  issue  will  be  considered  in  the  next  section,  in  the  context  of 
effective  parameter  extraction. 

Note  that  in  this  geometry  and  for  TM  incidence,  the  magnetic  field  is  parallel  to  the  fiber 
and  thus  has  a  nonzero  magnetic  flux  through  the  resonator  loop.  By  contrast,  arrays  of  thin 
SRRs  are  usually  printed  onto  a  substrate  so  that  the  resonator  loops  are  in  the  same  plane  as 
the  array;  experimental  characterization  of  one  to  a  few  layers  of  arrays  in  these  cases  can 
only  use  light  normally  incident  to  the  substrate  as  such  samples  are  too  thin  and  would  be 
strongly  diffracting  in  side  incidence.  The  magnetic  field  is  then  parallel  to  the  plane  of 
the  resonator  loops,  has  zero  magnetic  flux  and  is  unable  to  excite  a  magnetic  response. 
Instead,  the  electric  field  is  used  to  generate  currents  which  couple  to  the  magnetic  resonance 
[17]  while  transmission  features  due  to  a  magnetic  response  can  be  seen  in  these  experiments 
they  do  not  correspond  to  a  negative  permeability  [18]. 


3.  Characterization 


The  transmittance  of  our  samples  between  0.1-1  THz  was  measured  experimentally  via  THz 
time  domain  spectroscopy  [19,  20].  The  transmittance  results  for  the  TM  and  TE  waves  are 
shown  in  Figures  2(a)-(b)  for  the  two  fiber  orientations.  Under  TM  illumination,  our  arrays 
exhibit  distinct  magnetic  resonances  between  0.3 -0.4  THz,  depending  on  fiber  orientation. 
Under  TE  illumination,  they  possess  high-pass  filtering  behavior  (without  magnetic 
resonances),  analogous  to  the  response  of  metallic  wire-grid  structures  [21]. 

We  modeled  our  geometries  with  the  finite  element  solver  COMSOL,  using  a  periodic  2- 
dimensional  array  of  125  pm  Zeonex  squares,  coated  with  a  250  nm  silver  layer  on  three  sides 
(using  a  Drude  model  [15]  for  s)  and  unit  cell  size  350  pm.  A  polynomial  fit  to  experimental 
data  for  the  refractive  index  of  Zeonex  was  used  in  our  numerical  simulations.  There  was 
excellent  agreement  with  experiment  in  both  fiber  orientations  and  polarizations  [Figures  2(a)- 
(b)].  To  demonstrate  the  observed  resonances  are  indeed  magnetic,  we  consider  in  more  detail 
the  fields  inside  the  resonator:  Figure  2(d)  shows  the  relative  phase  of  the  simulated  magnetic 
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Fig.  4.  Real  and  imaginary  parts  of  the  retrieved  (a)  p,  (b)  8,  and  (c) 
magnetoelectric  coupling  coefficient  ^  for  the  symmetric  and  asymmetric  fiber 
orientations,  under  TM  polarization.  Note  that  for  symmetric  resonator 
orientations  ^  is  zero. 

field  at  the  center  of  the  resonator  with  respect  to  the  incident  field  at  that  point,  indicating  a 
near-7t  phase  shift  centered  around  resonance  at  0.36  THz.  Figure  2(e)  shows  a  color  plot  of 
the  simulated  magnetic  field  incident  under  TM  polarization  for  the  symmetric  resonator  (i) 
far  below  resonance  at  0.15  THz  (no  magnetic  response),  (ii)  just  below  resonance  at  0.34 
THz  (in-phase  magnetic  response)  and  (iii)  just  above  resonance  at  0.37  THz  (out-of-phase 
magnetic  response). 

We  extracted  p  and  s  of  the  metamaterial  arrays  from  simulations  under  TM  polarization 
using  the  retrieval  procedure  presented  in  Ref.  [16],  taking  into  account  the  bianisotropic 
behavior  of  our  asymmetric  fiber  array.  The  complex  scattering  matrix  elements  which 
describe  the  reflection  (Sn  and  S22)  and  transmission  (Si2  and  S2i)  of  the  field  incident  along 
the  x-axis  from  either  side  of  the  array  were  obtained  from  the  finite  element  calculations.  In 
the  symmetric  case  there  is  no  magnetoelectric  coupling,  Sn=S22  and  Si2=S2i;  in  the 
asymmetric  case  however  magnetoelectric  coupling  leads  to  Sn^S22  and  Si2=S2i.  Knowledge 
of  the  scattering  matrix  coefficients  and  metamaterial  thickness  enable  the  extraction  of  the 
refractive  index  n ,  from  which  c  and  p  and  the  magnetoelectric  coupling  coefficient  £,  can  be 
obtained.  The  results  are  presented  in  Figures  3(a)-(c).  The  symmetric  resonator  possesses  a 
minimum  p  =  -2+2.2/  with  no  magnetoelectric  coupling,  whereas  the  asymmetric  resonator 
displays  a  shallower,  red-shifted  magnetic  response  due  to  magnetoelectric  coupling,  with  a 
minimum  p  =  -0.9+ 1.2/.  Correspondingly,  s  exhibits  antiresonant  behavior  in  the  former  case, 
and  resonant  behavior  in  the  latter  [16,  22]. 


Note  that  in  this  work  we  have  measured  the  transmittance  and  extracted  the  effective 
parameters  of  a  single  metamaterial  layer  composed  of  periodic  sub -wavelength  resonators. 
Simulations  of  single  slabs  of  homogenized  material  with  the  extracted  s,  p  and  £,  indeed 
reproduce  the  same  transmission  and  reflection  curves  as  simulations  of  the  full  system  (not 
shown).  It  is  known  [23]  that  when  considering  stacks  of  multiple  layers,  these  parameters 
differ  from  the  single  layer  case  due  to  coupling  effects  between  adjacent  layers.  It  follows 
that  the  effective  s  and  p  depend  on  whether  the  layer  lies  inside  or  on  the  edge  of  an  arbitrary 
stack,  and  effective  parameters  of  stacked  layers  will  also  depend  on  the  way  layers  are 
stacked  (e.g.  distance  between  layers).  In  such  a  multiple-layer  case,  s  and  p  can  be  obtained 
by  analyzing  scattering  matrices  of  two-  and  three-  layer  systems  [23].  Such  an  analysis 
becomes  relevant  in  the  case  of  a  multi-layer  system,  and  is  beyond  the  scope  of  this  paper. 
Furthermore,  it  is  important  to  note  that  parameters  extracted  using  the  above  method  come 
from  matching  of  boundary  conditions,  but  do  not  necessarily  reproduce  the  intrinsic 
electromagnetic  response  within  the  metamaterial  [24]. 

4.  Conclusion 

In  conclusion,  we  have  presented  a  procedure  for  making  THz  metamaterials  with  a  negative 
effective  p  near  resonance.  The  excellent  agreement  between  experiment  and  simulations 
demonstrates  the  very  high  quality  of  the  fabricated  structures.  The  resonant  frequency  can  be 
tuned  simply  by  scaling  the  structure,  and  can  be  extended  to  higher  frequencies.  Drawing 
from  a  preform  allows  the  fabrication  of  fibers  of  -100  pm  diameter,  however  smaller 
diameters  (-10  pm)  can  be  obtained  with  standard  fibers  spinning  techniques  [25],  which 
would  move  the  resonant  frequency  up  to  -10  THz.  Further  optimization  may  allow  a 
reduction  in  size  down  to  1  pm,  allowing  the  fabrication  of  resonators  into  the  mid-IR. 
Whereas  for  THz  resonators  the  metal  thickness  is  far  smaller  than  the  wavelength  and  surface 
roughness  is  not  an  issue,  more  care  would  need  to  be  taken  at  shorter  wavelengths.  While  in 
our  fabrication  setup  we  use  a  batch  process  with  a  closed  sputtering  chamber,  in-line  coating 
with  a  continuous  operation  from  preform  to  metamaterial  fiber  is  possible,  with  the  unique 
advantage  of  being  able  to  fabricate  metamaterials  with  a  continuous  process  in  very  large 
quantities,  literally  by  the  kilometer.  Combined  with  the  previously  presented  fabrication  of 
fiber-based  negative  8  materials  [8]  this  will  enable  the  development  of  woven  negative  index 
materials,  as  well  as  the  fabrication  of  magnetic  surface  plasmon  waveguides  and 
sub  wavelength  waveguides  [26].  Advanced  weaving  techniques  will  allow  the  production  of 
materials  with  gradients  of  both  positive  and  negative  s  and  p,  providing  the  opportunity  to 
demonstrate  three  dimensional  cloaks,  as  well  as  hyperlensing  from  the  THz  to  the  mid-IR. 
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